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ABSTRACT 

We study the evolution of stellar mass in galaxies as a function of host halo mass, using the "MPA" 
and "Durham" semi-analytic models, implemented on the Millennium Run simulation. The results 
from both models are similar. We find that about 45% of the stellar mass in central galaxies in 
present-day halos less massive than ~ 10 12 /i _1 M Q is already in place at z ~ 1. This fraction increases 
to ~ 65% for more massive halos. The peak of star formation efficiency shifts toward lower mass halos 
from z ~ 1 to z ~ 0. The stellar mass in low- mass halos grows mostly by star formation since z ~ 1, 
while in high-mass halos most of the stellar mass is assembled by mergers. These trends are clear 
indications of "halo downsizing" . We compare our findings to the results of the phenomenological 
method developed by Zheng, Coil & Zehavi (2007). The theoretical predictions are in qualitative 
agreement with these results, however there are large discrepancies. The most significant one concerns 
the amount of stars already in place in the progenitor galaxies at z ~ 1, which is about a factor of 
two larger in both semi-analytic models. We also use the semi-analytic catalogs to test different 
assumptions made in that work, and illustrate the importance of smooth accretion of dark matter 
when estimating the mergers contribution. We demonstrate that methods studying galaxy evolution 
from the galaxy-halo connection are powerful in constraining theoretical models and can guide future 
efforts of modeling galaxy evolution. Conversely, semi-analytic models serve an important role in 
improving such methods. 

Subject headings: cosmology: galaxies - cosmology: theory - galaxies: statistics - galaxies: evolution 
— galaxies: halos - large-scale structure of universe 



1. INTRODUCTION 

In the current paradigm of structure formation, galax- 
ies form within cold dark matter halos. The formation 
and evolution of these halos are dominated by gravity 
and can be well predicted from high-resolution cosmo- 
logical numerical simulations and analytic models. The 
assembly of the stellar content of galaxies, however, is 
governed by more complex physics, and the relation be- 
tween galaxies and dark matter halos and the detailed 
physical processes of galaxy formation and evolution are 
only partially understood. 

A useful approach to explore galaxy formation within 
dark matter halos is the Se mi-Analytic Mo deling 
(SAM) of gal a xy formation Ce.g.. ICole et al 111994 l2fM 
Benson et all 120031: iCroton et al.l 120061: iBower et al.l 
20061 ). In such models, halos identified from high resolu- 
tion iV-body simulations are "populated" with galaxies 
using analytical prescriptions for the baryonic evolution. 
Within the SAM approach, galaxies change as the orig- 
inal stars evolve and new stars form. They also change 
their stellar content and increase their mass by merging 
with other galaxies. Different feedback or pre-heating 
mechanisms, such as those caused by star formation, 
active galactic nuclei, or the photo-ionizing ultra-violet 
background, also impact at different stages of a galaxy's 
life and are implemented in the models at different lev- 
els. These models have been successful in reproducing 
several measured properties including the galaxy lumi- 
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nosity and stellar mass functions (see e.g. .ICroton et al.l 
2006; B ower et alJl2006t iFontanot et al.ll2009t ). 

Different phenomenological methods have been 
developed to connect galaxies with dark mat- 
ter halos. One commonly used approach is the 
Halo Occupation Distrib ut ion fr amework (HOD, 
e. g. iJing. Mo. fc Boernerl 119981: iPeacock fc Smithl 



200a ISeliakl 120001: iScoccim arro et al.l [200lT 



Bcrlind & Weinbcr id 120021 ) . which characterizes the 



relationship between galaxies and halos in terms of the 
probability distribution, P(N\M), that a halo of virial 
mass M contains iV galaxies of a given type, together 
with the spatial and velocity distributions of galaxies 
inside halos. The HOD parameters are constrained 
using galaxy clustering measurements from large galaxy 
surveys and theoretically known halo clustering. Similar 
approache s include the Conditional Luminosi ty Function 
(CLF, see lYang. Mo. fc van den Boschll200l . which de- 
scribes the average number of galaxies as a function of lu- 
minosity that reside in a halo of mass M, and abundance 
matching schemes dConroy. Wechsler fc Kraytsovl 120061: 
Conrov fc Wechsler Il2009t iMqster et al.ll2010HGuo et al l 
2010t iNeistein et al.l 1201 Hi iRodriguez-Puebla et al.l 



20111 ). which monotonically connect galaxy luminosity 



or stellar mass to halo mass by matching the abundances 
of halos and galaxies. 

HOD models have been mostly used to learn about 
the relationship between galaxies and halos at a 
fixed epoch (e.g.. iBullock Wechsler. fc Somervillell2002t 
IZehavi et all 120051 120111: IZheng et al.l 120081 and ref- 
erences therein). Recent studies have started using 
them to also explore galaxy evolution by combining 
the inferred galaxy-halo connection at different red- 
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shifts with the evolution of dark matter halos pro- 
vided by theory dZheng et al I [20071: iWhite et all 120071: 



ISeo. Eisenstein fc Zehavif!2008F iWake et al ll2008l l20TTh . 

In particular, Zheng, Coil & Zehavi (2007; hereafter 
ZCZ07) develop a phenomenological approach to extract 
information about galaxy evolution from z ~ 1 to z ~ 
by performing HOD modeling of two-point correlation 
functions of DEEP2 and SDSS galaxies. With the in- 
ferred galaxy-halo connection at two redshifts, they es- 
tablish an evolutionary link using the typical growth of 
dark matter halos obtained from numerical simulations. 

Even with the progress made in establishing the evo- 
lutionary link between galaxies and halos, our under- 
standing of the specifics of stellar mass growth within 
the dark matter halos is still far from complete. Galax- 
ies can grow their stellar mass by star formation, accre- 
tion of smaller satellite galaxies or major merging. It is 
important to quantify the contribution of all these pro- 
cesses in order to have a complete picture of the assembly 
history of galaxies within their host dark matter halos. 
ZCZ07 derive the mean stellar masses of central galax- 
ies at z ~ 1 and z ~ as a function of the present-day 
host halo mass. After roughly accounting for the contri- 
bution of merging of central and satellite galaxies, they 
infer the star-formation contribution to the stellar mass 
assembly. They find that in central galaxies located in 
relatively low-mass halos (~ 5 x l0 11 h~ 1 M ( ? ) ) the bulk 
of the stellar mass results from star formation between 
z ~ 1 and z ~ 0, while only a small fraction of stars 
formed since z ~ 1 in central galaxies of halos as massive 
as ~ 5 x 10 12 /i~ 1 M Q . For these massive halos, merg- 
ing becomes more important and constitutes the domi- 
nant contribution to the stellar mass growth (see Fig. 9 
in ZCZ07 for details). The results reflect the so-called 
"downsizing" star formation pattern in which the sites 
of active star formation shift from high-mass galaxies at 
early times to low er-mass systems at later epochs (e.g., 
iCowie et al.lll99 6). manifested in terms of halo mass. 

In this paper, we study the theoretical predictions for 
stellar mass evolution as a function of dark matter halo 
mass using SAM catalogs. One of the main advantages 
of the SAMs is that we can trace the full evolution of 
the individual galaxies within their dark matter halos, 
allowing an explicit study of the different processes that 
contribute to the build up of the stellar content of galax- 
ies. We compare and contrast these predictions with the 
results of ZCZ07. We gauge the potential of the phe- 
nomenological methods to constrain galaxy formation 
models, as well as test some of the assumptions of such 
methods. In particular, we check the validity of the sim- 
ple evolutionary approach presented by ZCZ07. 

For th ese purposes, we u s e two SAM catalogs, the 
"MPA" (ICroton et all [2001 IDe Lucia fc Blaizoti 120071) 
and "Durham" ([Bower et al.ll2006t iFont et aLll2008D cat- 
alogs produced fro m the Millennium R un cosmological 
N-hoAy simulation (Sp ringel et al.|[2005[) . We mainly fo- 
cus on the stellar mass growth as a function of halo mass 
since z ~ 1 for ease of comparison with ZCZ07. Nonethe- 
less, we briefly study also the stellar mass growth since 
z ~ 2 in the MPA catalog to explore the processes in- 
volved in galaxy formation at higher redshifts. These 
SAM catalogs have been previously use d to study the 
stellar mass evolution in galaxies (e.g.. |Guo fc White! 
120081: iStringer et all 120081: iFontanot et al J I2009D . How- 



ever, those studies focus on the integrated stellar mass 
to make a direct comparison to observations and did not 
investigate in detail its evolution as a function of halo 
mass. Most physical processes involved in galaxy forma- 
tion models depend explicitly on the mass of the dark 
matter halo. Additionally, the modeling of other ob- 
servables such as the galaxy-galaxy merger rate strongly 
rely on the pr ecise knowledge of th e galaxy-halo connec- 
tion (see e.g.. lHopkins et alll2010D . Thus it is physically 
meaningful and informative to study galaxy evolution as 
a function of halo mass. 

The paper is organized as follows: §2 describe the SAM 
catalogs that we use and the sample selection. In §3 we 
present and discuss our results for the growth of stellar 
mass as a function of halo mass. In §4 we compare our 
results to those found by ZCZ07 and we conclude in §5. 

2. MOCK CATALOGS AND GALAXY FORMATION 
MODELS 

In this work, we use the publicly av ailable mock galaxy 
catalogs produced with th e "MPA" (ICroton et alj|2006t 
| De Lucia fc Blaizot 2007) and "Durham" ( Bower et al.l 
l2006f ) semi-analytic models of galaxy formation, both 
based upon dark matter halo evolution in the Millennium 
Run simulation ([Springel et al.l 120051) . The Millennium 
Run followed the evolution of ~ 10 10 dark matter parti- 
cles in a ACDM Universe. The simulation uses a peri- 
odic box of 500/i~ 1 Mpc on a side with mass resolution 
per particle of 8.6 x 10 s h~ 1 MQ. The initial conditions of 
the simulation were generated with cosmological parame- 
ters obtained from the combined analysis of 2dFGRS and 
WMAP1 CMB data. The halos in the simulation were 
identified in each time step using a friend-of-friends algo- 
rithm with a linking length of 0.2 times th e mean particle 
separ ation. More details can be found in ISpringel et all 
[20051 

In the SAMs, galaxies are assumed to form at the cen- 
ter of dark matter halos. The evolution of the baryonic 
component of galaxies is modeled using simple but phys- 
ically motivated analytic prescriptions. These include 
radiative cooling of hot gas, star formation in the cold 
disk, supernova feedback, black hole growth and AGN 
feedback through the "quasar" and "radio" epochs of 
AGN evolution, metal enrichment of the inter-galactic 
and intra-cluster medium, as well as galaxy morphology 
shaped through mergers and merger-induced starbursts. 

As mentioned above, these models are aimed at re- 
producing integrated galaxy observables. To that ef- 
fect, the SAMs recover reasonably well t he galaxy lu- 
minosity function in diff erent bands (e.g., ICroton et al.l 
[2001: IB ower et al"1l2006f) as well as th e stellar mass func- 
tions for a range of redshifts (e.g., IBower et al.l | 2006h 
Fontanot efaLl 120091: but see also ILi fc White! 12009; 
Marchesini et al.l [2009) . The predictions for the frac- 
tion of blue and re d central galaxies , however, are 
not fully correct (see lBaldrv et al.l I2006T ) . a nd modifica- 
tions to the treatment of gas cooli ng (e.g.. iViola et all 
[2001) and the AGN feedback (e.g., IBaldrv et all 120061) 
have been proposed to alleviate these discrepancies. 



red satellites dWeinmann et all 120061: 


Coil et all 120081: 


Kang & van den Bosch 2008: Do Lucial 


20091: iKim et all 


2009). Note that satellite galaxies in these models are 



treated differently than central galaxies. Only central 
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galaxies accrete new material by cooling from the hot 
atmosphere of their halo, direct infall of cold gas, and 
the merging of satellites. Since no new material accretes 
onto satellites, their star formation ends when the cold 
gas is exhausted (see lCroton et all 12006'). Satellites are 
also affected by different environmental processes that 
change their properties, which ar e not fully implem ented 
in these models. For example, iKim et all ()2009f ) sug- 
gest that satellite-satellite mergers and tidal dissolution 
of satellites need to be included to better match the mea- 
sured luminosity dependence of galaxy clustering. 

Even though different SAMs adopt similar analytical 
prescriptions to treat the physical processes involved in 
galaxy formation and evolution, there are significant dif- 
ferences in the way the MPA and Durham SAMs deal 
with specific processes, such as the cooling of gas, the cut- 
off bl ack hole mass for AGN feedback (see Stringe r et al.1 
2008) and the dynamical treatm ent of the "orphan" 
satellite galaxies (fGao et al.ll2004fl . Also, the two models 
are based on different halo merger trees which are con- 
stru cted in the post-processing stage of th e simulation 
(see lBower et alj|2006t fPe Lucia fc Blaizotl 120071) . Even 
though the halo merger trees are checked to be statis- 
tically compatible (De Lucia, private communication), 
differences could arise at the level of individual galaxies, 
affecting some galaxy properties. 

The full SAM galaxy catalogs used in this work0 con- 
tain information for about eight million galaxies brighter 
than M r = — 17. Several properties are available for each 
of these galaxies, including positions and velocities, mag- 
nitudes in several band passes (Johnson, Busher, 2MASS 
as well as the 5 SDSS bands), stellar mass, and mass of 
the dark matter halo in which the galaxy is located. It 
is important to note that these two SAMs assume differ- 
ent initial mass functions of stars . The MPA mo del as- 
sumes a Chabrier mass function (|Chabrierl l2 00l. while 
the Durham model uses the Kennicutt one ()Kennicuttl 
1983). In order to make a direct comparison between 
the models, we consistently t ransform both t o a "diet" 
Salpeter initial mass function (jBell et al.ll2003| ) , used also 
in ZCZ07. 

We select from each catalog a random sample of 
250,000 present-day central galaxies. For compatibility 
with the ZCZ07 results, we use in fact the z ~ 0.1 snap- 
shot of the catalogs, as this is approximately the value 
of the median redshift for SDSS galaxies. We hereafter, 
however, loosely refer to it as z ~ 0, in comparison to the 
evolution since z ~ 1 and z ~ 2. For each present-day 
galaxy, we identify all its progenitors at z ~ 1 (and also 
at z ~ 2 in the MPA SAM) by following the merger tree 
of each present-day central galaxy. These progenitors can 
be central galaxies in dark matter halos, satellites located 
in subhalos or "orphan" galaxies, i.e. satellite galaxies 
whose parent dark matter subhalo was destroyed below 
the resolution limit of the simulation by ti dal stripping 
and truncation (see, e.g., iGao et all 120041 ) . We define 
the main progenitor of a z ~ galaxy as the central 
galaxy located in the most massive dark matter halo at 
the corresponding higher redshift {z ~ 1 for both cat- 
alogs and z ~ 2 for the MPA only). We have verified 
that defining the main progenitor as the most massive of 

4 The S AM galaxy catalogs can be downloaded from 
http: / /www. g-vo.org/Millennium/ 



the m erger tree main branch (e.g., iDe Lucia fc Blaizotl 
l2007f ) does not change our results. Once all progenitors 
are identified, we can quantify the different contributions 
to the stellar mass of central galaxies, namely the contri- 
butions from smaller central galaxies and from satellites 
that merge into the central galaxies. The present-day 
stellar mass in central galaxies that does not come from 
the different merger processes arises then from new star 
formation during that time period. 

3. STELLAR MASS GROWTH AS A FUNCTION OF HALO 

MASS 

In this section we first evaluate the star formation ef- 
ficiency of galaxies as a function of host halo mass. We 
then study the evolution of stellar mass in central galax- 
ies in the semi-analytic models. We also analyze the dif- 
ferent contributions to the stellar mass growth of central 
galaxies since z ~ 1. Finally, we briefly investigate the 
stellar mass growth since z ~ 2. 

3.1. Star Formation Efficiency 

A first fundamental quantity to investigate is the star 
formation efficiency (SFE) as a function of halo mass. 
Here, we define the SFE as the total stellar mass in 
central galaxies (M sta r) divided by the baryon mass 
Mb = fbM, where fb is the global baryon fraction and M 
is the mass of the dark matter halo. This SFE represents 
the integrated value from the redshift of formation to the 
epoch in consideration, reflecting the fraction of baryons 
associated with halos that are converted into stars in that 
time period. 

Figure [T] shows the SFE at z ~ 1 and z ~ as 
a function of halo mass predicted by the MPA SAM. 
The baryon fraction adopted in the MPA SAM is fb — 
0.17, though the baryon fraction in any individual halo 
may vary from the global value. The SFE at both 
epochs has a similar shape, peaking at some charac- 
teristic mass and dropping towards both low and high 
mass ends. At z ~ 0, the SFE peaks for halos of mass 
~5x 10 11 /i _1 M©, while at z ~ 1 the peak corresponds 
to more massive halos (~ 10 12 /j Mq). This trend 
is an indica t ion o f the "downsiz i ng" phenomena (e.g., 
Cowie et al.l 119961: Uuneau et al.l | 2005t iFontanot et al.l 
20091 : lAvila- Reese &: Firmanil 1201 ll ). shown specifically 
here as a function of halo mass. The SFE at the peak is 
~ 23% and ~ 18% at z ~ and z ~ 1, respectively. The 
SFE results for the Durham model are similar. 

The drop of SFE at the low- mass end is associated with 
the availability of cold gas in these halos, which can be 
affected by photo-ionization heating and star formation 
feedback. The drop of SFE above the characteristic halo 
mass is related to gas accretion becoming less efficient 
due to the high virial temperature and AGN feedback. 
Additionally, we are only considering here central galax- 
ies, and in high-mass halos the stellar mass contributions 
from satellite galaxies can be substantial. 

3.2. Stellar Mass Growth of Central Galaxies 

Figure [2] shows the mean stellar mass in central galax- 
ies at z ~ (thick lines) and that of their z ~ 1 main 
progenitors (thin lines) as a function of the present-day 
halo mass. The results for the MPA catalog (solid lines) 
and the Durham catalog (dashed lines) are similar. At 
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Fig. 1. — Star formation efficiency at 2 ~ (solid line) and at 
2 ~ f (dashed line) as a function of halo mass for the MPA SAM. 

both redshifts, the stellar mass of the central galaxy in- 
creases rapidly with halo mass at the low mass end and 
slowly at the high mass end. The transition halo mass is 
approximately 10 12 /i _1 M© at z ~ and 2 x 1O 12 /i _1 M 
at z ~ 1. 

The bottom panel shows the ratio between the stellar 
mass at z ~ 1 and z ~ 0, representing the fraction of 
the z ~ central galaxy stellar mass already in place in 
the progenitor central galaxies at z ~ 1, as a function of 
present-day halo mass. For example, for a present-day 
halo of mass 5 x 10 11 /i- 1 M Q , on average - 50% (~ 40%) 
of the stellar mass in the central galaxy is already in 
place in the z ~ 1 main progenitor central galaxy for the 
MPA (Durham) model. The ratio gradually increases 
to ~ 65% for halos with mass ~ a few x 10 12 h~ 1 M.Q, 
and it starts decreasing towards the highest halo masses 
probed in this work. The decrease of stellar mass ratio 
at the highest halo masses is in accord with the predic- 
tions for hierarchical assembly of massive gala xies (e.g., 
iDe Lucia et al.ll2006t iDe Lucia & Blaizotll2007t ). 

3.3. Different Contributions to the Stellar Mass Growth 

The stellar mass in galaxies grows as a consequence of 
internal star formation or external infall of material (ma- 
jor and minor mergers). For the former, observatio nal 
estimates fe.g.. lSalim et~aT1[2007t iNoeske et al.ll2007h of- 
ten come from measuring the average star formation rate 
as a function of stellar mass and time. The amount of 
stellar mass gained through accretion is more difficult to 
measure directly. It is often estimated by simply taking 
the difference between the growth due to star formation 
and the total stellar mass at present. 

In a SAM based on an A-body cosmological simula- 
tion, we have the full merger history of dark matter halos 
and galaxies. Thus, it is possible to track the complete 
evolution of the stellar mass due to both mergers and 
star formation as a function of the halo mass. In par- 
ticular, following ZCZ07, we account for four different 
components of the assembly of stellar mass in central 
galaxies: stars in place in the progenitor central galax- 
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Fig. 2. — Mean stellar mass in z ~ central galaxies and their 
z ~ 1 progenitor central galaxies as a function of the present-day 
halo mass. Top panel: The mean stellar mass in central galaxies at 
2 ~ (thick lines) and that of their z ~ 1 main progenitors (thin 
lines) as a function of the present-day halo mass, predicted by the 
MPA (solid lines) and Durham (dashed lines) SAMs. Bottom panel: 
The ratio between the stellar mass at z ~ 1 and z ~ for the MPA 
(solid line) and Durham (dashed line) SAMs, which represents the 
fraction of z ~ central stellar mass that is in place in the z ~ 1 
progenitor central galaxies. 

ies at z ~ 1, stellar mass from smaller central galaxies 
that merge to the central galaxy, stellar mass from any 
satellites (of the main progenitor or other smaller central 
galaxies) that merge with the central galaxy, and recent 
star formation. 

Figure [3] presents the SAM predictions for these dif- 
ferent contributions to the z ~ central galaxy stellar 
mass as a function of the present-day halo mass, for both 
the MPA and Durham models. In each panel, the bot- 
tom curve (marked as "A" ) denotes the fraction of stellar 
mass in place in the z ~ 1 progenitor central galaxies and 
is essentially the same ratio shown in the bottom panel 
of Figure [2] The area between curves "A" and "B" in- 
dicates the contribution of smaller central galaxies that 
have merged into the main progenitor. The area between 
curves "B" and "C" shows the contribution from satellite 
galaxies in all progenitor halos at z ~ 1 that end up in 
the z ~ central galaxies. Consequently, the remaining 
stellar mass (from curve "C" to the top of the plot) arises 
from star formation since z ~ 1. 

In the SAM datasets, the amount of star formation 
since z ~ 1 can also be directly obtained by integrating 
the star formation rate as a function of time. We do 
a crude calculation of this in the MPA SAM for three 
specific present-day halo masses, as a sanity check, by 
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Fig. 3. — Different contributions to the present-day central galaxy stellar mass as a function of present-day halo mass shown for the MPA 
(left panel) and Durham (right panel) SAMs. The first curve from the bottom (denoted as "A") is the stellar mass already in place in 
the progenitor central galaxies at z ~ 1. The area between lines "A" and "B" indicates the contribution of the smaller central galaxies 
that merge with the progenitor central galaxies. The area between curve "B" and "C" denotes the contribution of satellite galaxies. The 
remainder (from curve "C" to the top of the plot) is the stellar mass gained through star formation since z ~ 1. The full circles in the left 
panel are the contribution of star formation in the main progenitors obtained directly by integrating the star-formation rate over time (see 
text for details). 



integrating the star-formation rate for a subsample of 
main progenitors over 10 snapshots since z ~ 1 (shown 
as filled circles in the left panel of Figure ■ This gives 
a sense of the total star formation since z ~ 1, since the 
contribution from smaller central progenitors is minor 
and the satellites do not contribute to the star formation 
in these models (see § 2). The first point appears to be a 
bit below the expected value (corresponding to a larger 
star formation contribution), reflecting the approximate 
nature of this calculation. 

The overall trends of stellar mass growth from the two 
SAMs are similar. It is evident that central galaxies in 
small halos grow mostly by star formation since z ~ 1, 
while the star formation contribution is small in high 
mass halos. This could be explained by the fact that 
small halos are almost completely assembled by z ~ 1 , so 
the new stellar mass must come from star formation. For 
intermediate halo masses, the contribution from merging 
becomes more significant, but star formation still con- 
tributes most of the stellar mass since z ~ 1. For high- 
mass halos, the contribution from mergers dominates 
and star formation is negligible. This overall behavior 
is another manifestation of the "downsizing" effect, in 
this case referring to the fact that more massive galaxies 
form the bulk of their stars earlier tha n smaller galax- 
ies ( " anti-hierarchical" ass embly; seee.g.jDe Lucia et al.l 
2006; iStringer et al]l200l iFontanot et al.ll2009t ). 

Some differences between the two SAM models are ap- 
parent. The stellar mass already in place in the z *~ 1 
progenitors was compared in the bottom panel of Fig- 
ure[2l With regard to the other components, the Durham 
model predicts a larger contribution from smaller central 
progenitors in halos with present-day mass larger than 
~ 10 12 /i _1 M Q , while the MPA model produces a slightly 



larger contribution from satellite galaxies in low-mass 
halos. These discrepancies likely reflect the differences 
in the galaxy formation prescription in the two models 
as well as the differences in the underlying halo merger 
trees. Another contributing factor is the different treat- 
ment of the dynamics of "orphan" galaxies (De Lucia, 
private communication) . 

3.4. Stellar Mass Growth Since z ~ 2 

In this section we explore the SAM predictions for 
stellar mass growth from higher redshifts, to gain fur- 
ther insight on the processes contributing to galaxy 
evolution and serve as a reference point for model- 
ing higher-redshift observations. The interpretation of 
these predictions is complex since observations currently 
do not provide consistent indicators of stellar masses 
and star formation rates at these redshifts (see e.g., 
iConrov &; Wechslerl 120091 §3.5 for a discussion). Also 
clustering data at z ~ 2 is relatively scarce for ZCZ07- 
like analyses (but see I Wake et al.|[20lTl for a first attempt 
along these lines). 

We study the stellar mass growth as a function of halo 
mass since z ~ 2 in the MPA catalog, using a smaller 
sample of 160, 000 galaxies. The SFE at z ~ 2 as a 
function of halo mass exhibits a similar shape and peak 
location as the z ~ 1 results (shown in Fig. |TJ) . The over- 
all amplitude is slightly smaller, with a SFE of ~ 16% at 
the peak of the distribution. Our results for the growth 
of stellar mass are presented in Figured] The top panel 
shows the mean stellar mass in central galaxies at z ~ 
(thick line) and that of their z ~ 2 main progenitors 
(thin line), as a function of the present-day halo mass. 
The bottom panel presents the different contributions to 
the present-day central galaxy stellar mass. The differ- 
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Fig. 4. — Top panel: Mean stellar masses in central galaxies at 
z ~ (thick line) and that of their z ~ 2 main progenitors (thin 
line), as a function of the present-day halo mass, predicted by 
the MPA SAM. Bottom panel: The different contributions to the 
present-day central galaxy stellar mass. Curves are the same as in 
Fig. [3] except now referring to stellar mass evolution since z ~ 2. 

ent curves in this plot are the same as those in Figure 
but referring now to z ~ 2, with curve "A" representing 
the stars already in place in the central galaxies at z ~ 2. 

The trends are qualitatively very similar to the ones 
seen for the z ~ 1 progenitors in the previous sections. 
As expected, there is significantly less stellar mass al- 
ready in place in the main progenitor at z ~ 2. The 
contribution from the smaller central galaxies is more 
significant than that for the z ~ 1 case, while the con- 
tribution of the satellites is roughly the same. As less 
stellar mass is in place at z ~ 2, the contribution of star 
formation from z ~ 2 to the present-day is substantial, 
at all halo masses; it is about 10% for central galaxies 
in the highest mass halos probed and higher for those in 
lower mass halos. Note that the "downsizing" pattern is 
also evident here. 

4. COMPARISON WITH PHENOMENOLOGICAL 
APPROACH 

Zheng, Coil, & Zehavi (2007; ZCZ07) perform HOD 
modeling of the luminosity-dependent projected two- 



point correlation function for DEEP2 and SDSS galaxies, 
at z ~ 1 and z ~ 0, respectively. They infer the relation- 
ship between central galaxy luminosity and halo mass at 
these two redshifts and establish an evolutionary link by 
using the typical growth of dark matter halos obtained 
from numerical simulations. Stellar masses are derived 
from the galaxies luminosity and color. As a proof of con- 
cept, they estimate the evolution of galaxy stellar mass 
as a function of host halo mass. An approximate method 
is used to estimate the different contributions of merg- 
ers and star formation to the growth of central galaxies 
stellar mass. In this section, we use the SAM results to 
gauge the potential of such phenomenological methods 
to constrain galaxy formation and evolution models. We 
first compare the stellar mass growth in central galax- 
ies inferred by ZCZ07 from DEEP2 and SDSS galaxy 
clustering with that predicted by the SAM models. We 
then examine the validity of the assumptions used in the 
ZCZ07 approach. 

4.1. Star Formation Efficiency and Galaxy Stellar Mass 

Growth 

We first examine the star formation efficiency and its 
evolution from rcdshift z ~ 1 to obtained using the 
MPA SAM model (our Fig. [TJ and in ZCZ07 (their 
Fig. 10). Both approaches produce similar general trends 
with halo mass, exhibiting peaked distributions at the 
two redshifts, and at similar halo masses. The ZCZ07 
results also exhibit the halo "downsizing" effect, with the 
SFE peak shifting to a higher mass at the higher redshift. 
The MPA SAM and ZCZ07 also find comparable values 
for the maximum SFE at z ~ 0. However, at 2 ~ 1, 
the SAM shows an overall higher SFE than computed in 
ZCZ07 (peaking at 18% and 12%, respectively), which 
translates into a larger amount of stars by that redshift. 

Figure[5]shows the SAM predictions (solid lines) for the 
stellar mass as a function of halo mass and the results ob- 
tained by ZCZ07 (dashed lines), over the halo mass range 
they probe. Although the SAMs and ZCZ07 methods 
produce similar trends, there are important quantitative 
disagreements between them. The main difference is that 
the SAMs predict many more stars already in place at 
z r~j 1 in the progenitor central galaxies compared to the 
ZCZ07 results. The differences are especially pronounced 
at medium and low halo masses. This may be related 
to the known fact that the SAM s produce too many 
M+ g alaxies at high redshift (e.g., IKitzbichler fc White! 
I2007f l . At z ~ 0, on the other hand, the agreement is 
quite good for low mass halos, while for high-mass halos 
(larger than ~ 1O 12 /i _:l M ) the SAMs seem to underpre- 
dict the stellar mass in central galaxies. The AGN "radio 
mode" feedback be comes important on these mass scales 
(|Bower et al.ll2006f ). This suggests that the SAMs might 
be overestimating the strength of the feedback. 

As for the fraction of stellar mass in place in the z <~ l 
progenitor central galaxies, the SAMs prediction is about 
twice that inferred by ZCZ07, as shown in the bottom 
panels of Figure El The error bars in the SAMs denote 
the la scatter around the mean. Note that this discrep- 
ancy could be less dramatic, as there might be ~ 25% 
underestimation in the ZCZ07 calculation of the stellar 
mass at z ~ l due to DEEP2 red galaxy incompleteness 
(see ZCZ07 for details). This effect is shown by the dot- 
ted lines in the bottom panels. Even with this potential 
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Fig. 5. — Top panels: Comparison of the MPA (left, solid lines) and Durham (right, solid lines) SAMs with the ZCZ07 results (dashed 
lines) for the mean stellar mass in central galaxies at z ~ and z ~ 1, as a function of present-day halo mass. Bottom panels: Comparison 
of the SAMs (solid lines) and ZCZ07 results for the ratio of central galaxy stellar mass at 2 ~ 1 and z ~ as a function of present-day 
halo mass. Error bars in the SAMs predictions are the lcr scatter around the mean. The dashed line (denoted as "ZCZ07") represents the 
standard ZCZ07 result. The dotted line (labeled as "ZCZ07+") is the conservative estimation, assuming a 25% correction of the stellar 
mass in the DEEP2 sample (see text for details). 



correction, the discrepancy is significant. 

We note that these discrepancies are present at roughly 
the same level for both SAMs, indicating that their cause 
is of a more fundamental origin not reflected in the differ- 
ences between the two models. It is also worth mention- 
ing that the major discrepancy between the SAMs and 
ZCZ07 predictions for the stellar mass growth appears 
to be present already at redshift 2. From the results pre- 
sented in §3.4, we see that the amount of stars in place 
in central galaxies by z ~ 2 predicted by the MPA SAM 
is of the order of the phenomenological results for the 
amount of stars in place by z ~ 1. 

The top panels of Figure [6] compare the predictions 
of the SAMs (solid lines) to the ZCZ07 results for the 
total stellar mass acquired through merging of smaller 
central and satellite galaxies on top of that already in 



place at z ~ 1, normalized by the final stellar mass at 
z ~ 0. For ZCZ07, we plot both the standard estimation 
(dashed line in each panel) and the conservative estimate 
including the possible 25% correction of the stellar mass 
at z ~ 1 (dotted line). Here the agreement is better, 
especially at high halo masses. The difference between 
these total contributions and the stellar mass at present- 
day in halos of a given mass is essentially the contribution 
arising from star formation since z ~ 1. 

Finally, the bottom panels of Figure [5] contrast the in- 
dividual contributions to the stellar mass assembly in the 
SAMs (solid lines) and ZCZ07 (dashed lines), showing 
the relative contributions to the final central galaxy from 
mergers of smaller central galaxies (thick lines) and satel- 
lite mergers (thin lines). The individual contributions 
from these models are strikingly different. In ZCZ07, the 




Fig. 6. — Top panels: Predictions of the SAMs compared to the results obtained by ZCZ07 for the total stellar mass in central galaxies 
that was acquired through merging, in addition to that already in place at z ~ 1. The stellar mass is normalized by the amount of stellar 
mass at z ~ and presented as a function of present-day halo mass. The left and right panels show the comparison to the MPA and 
the Durham SAMs, respectively (solid lines). The dashed lines are the standard ZCZ07 result and the dotted lines denote the effect of a 
tentative 25% correction of the stellar mass at Z ~ 1 (see text for details). Bottom panels: Comparison between the SAMs (solid lines) and 
ZCZ07 (dashed lines) for the individual contributions from mergers of smaller central galaxies (thick lines) and satellites (thin lines) to the 
stellar mass in the final central galaxy. 



main merger contribution to present-day central galax- 
ies comes from mergers of the smaller central progenitor 
galaxies, while in the SAMs a larger contribution comes 
from satellite galaxies. As mentioned already in §3.3, the 
contribution from smaller central progenitor galaxies in 
the MPA model is particularly minor, while it is some- 
what larger in the Durham model. The implications of 
these differences are discussed further below. 

4.2. Examining Assumptions in the ZCZ07 Approach 

In this section we use the SAM catalogs to test the 
validity of some of the assumptions adopted in the ZCZ07 
method. We note that these tests can be done despite 



the discrepancies found between the SAMs and ZCZ07. 

In ZCZ07, an evolutionary link is established between 
galaxy populations at the two epochs via the theoret- 
ically predicted growth of dark matter halos in which 
these galaxies reside. The first assumption we exam- 
ine is related to the statistical nature of any such rela- 
tion. ZCZ07 use the average relation between masses of 
present-day halos and that of their main progenitors at 
z ~ 1 , obtained by the PINOCCHIO code ((Monaco et al.l 
120021) . with no accounting for the intrinsic scatter. This 
could, in principle, impact the estimation of the stellar 
mass growth. 

Using the SAMs, where the full merger tree is known, it 



9 



o 



v 




H — I I I I I 



Full MPA SAM 
Averaged' MPA SAM 

1 1 1 — I I I I I 



_L 



1 

0.8 

.2 0.6 
*-> 

CO 

^ 0.4 
0.2 


10 18 1Q13 

M (h" 1 M ) [z~0] 

Fig. 7. — The MPA predictions for the mean stellar mass in 
central galaxies at z ~ and that of their z ~ 1 main progenitors 
as a function of the present-day halo mass, calculated using the full 
halo growth information (solid lines) and the "averaged" relation 
from the ZCZ07 approach (dashed lines). See text for details. 

is possible to test this. For each individual z ~ central 
galaxy, the SAM catalog provides its present-day stellar 
mass and host halo mass as well as its main progenitor 
stellar mass and host halo mass at z ~ 1. Based on these, 
we derive the the stellar mass of central galaxies and that 
of the progenitor central galaxies as a function of the 
present-day halo mass, which by construction includes 
the scatter in the halo mass and progenitor halo mass 
relation. With the SAM catalog, we can also follow the 
ZCZ07 procedure to obtain the above relations by using 
the average growth of dark matter halos. 

We perform such a test with the MPA catalog. Fig- 
ure [7] compares the stellar mass evolution as a function of 
present-day halo mass obtained with those two different 
procedures. The solid lines are the results obtained us- 
ing the individual halo growth information (presented al- 
ready by the solid lines in Fig. 12}, while the dashed lines 
denote the predictions using the average halo growth. 
There are only negligible differences between both re- 
sults, with a maximum deviation of about 5%, indicat- 
ing that the use of the average halo growth to connect 
galaxies at the two epochs is adequate. 

The other assumptions in ZCZ07 are related to the 
rough estimation of the merger contributions to the stel- 
lar mass growth of central galaxies, resulting in the 
dashed curves shown in the bottom panels of Figure [5J 
Computing these contributions in a statistical way is not 
straightforward. In order to estimate the stellar mass 
contribution from smaller central galaxies, ZCZ07 uti- 



lize, for a given dark matter halo mass, the fraction of 
halo mass already formed at z ~ 1 in its main progenitor. 
This halo has some fraction of stars already formed by 
that epoch. The assumption is that the ratio of central 
galaxy stellar mass to halo mass in the other progeni- 
tors is the same as that for the main progenitor and that 
they all will have merged into the main progenitor cen- 
tral galaxy by the present time. For example, a halo of 
mass 5 x lO 11 /! -1 ^/!© (at present) has assembled about 
70% of its mass by z ~ 1, while the MPA SAM predicts 
that the corresponding central galaxy has about 50% of 
the stars already in place by z ~ 1. Following the ZCZ07 
procedure, the remaining 30% of the halo should con- 
tribute 30/70 x 50% - 21% of the stars. Implicitly, this 
computation assumes a constant star formation efficiency 
(or that the merging halos are of comparable mass). Ob- 
taining the contribution of satellite galaxies to the final 
central galaxies is more difficult. ZCZ07 provide a crude 
estimation based on simplifying assumptions resulting in 
a linear extrapolation from zero at the low-mass end, 
where the satellite contribution is expected to be negli- 
gible, to 25% of the central galaxies contribution at the 
high-mass end, where the brightest satellite in each halo 
is assumed to have merged with the central galaxy (see 
their §6.3 for more details). 

We use the full halo and galaxy merging histories pro- 
vided in the SAMs to test these assumptions regarding 
the contribution from smaller central galaxies and from 
satellites. The top panels of FigureOshow for both SAMs 
the predicted contribution of smaller central galaxies to 
the final stellar mass obtained by applying the simple es- 
timation proposed by ZCZ07 (dashed lines). The exact 
contribution in the SAMs is shown by the solid lines. It 
is apparent that the approximate approach overestimates 
the "real" contribution by a large amount. Recall that 
the SAMs appear to predict a large excess of stars that 
already assembled in the central progenitors at z ~ 1 
(Fig. |SJ), which might translate to an overestimation of 
the estimated contribution. 

Another important issue is that the estimation assumes 
that all the remaining halo mass accreted carries with 
it the same fraction of stellar mass. This may be rea- 
sonable for the merging smaller halos, however, there 
is also a significant contribution from smooth "diffuse" 
accretion of dark matter particles to the final halo. In 
the Millennium Simulation, we obtain that this compo- 
nent accounts for ~ 30% of the final halo mass since 
z ~ 1 (see alsolGuo fc Whitel l2oM iFakhouri fc MallMfil : 
IGenel et al.ll2010l) . Even though smooth accretion carries 
baryons that will be added to the hot gas available for 
cooling ()De Lucia et al.l [2001) . it does not contribute to 
the existing stellar mass. When accounting for this, we 
find that the contribution from smaller central galaxies 
is considerably reduced (the dot-dashed lines in top pan- 
els of Fig. [5]) and is in better agreement with the SAMs 
predictions (particularly for the Durham model). This 
approximated contribution, however, still appears to be 
somewhat overestimated, especially in the MPA model, 
where the measured contribution of smaller centrals is 
tiny. 

As mentioned already, the simplified estimate also im- 
plicitly assumes a constant SFE for the merging halos, 
which in reality does vary with halo mass (see, e.g., Fig.[T] 
in this work and Fig. 10 in ZCZ07). This can affect the 
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Fig. 8. — Testing ZCZ07 assumptions regarding the contribution to the central galaxy stellar mass from mergers of smaller central 
galaxies and satellites since z ~ 1 with the MPA SAM (left) and Durham SAM (right). The contribution is normalized by the total amount 
of stellar mass at z ~ and presented as a function of present-day halo mass. Top panels: The solid line in each panel denotes the "exact" 
stellar mass contribution from the merging of smaller central galaxies. The predictions for this contribution using the ZCZ07 assumptions 
applied to the SAMs are shown as the dashed-lines (marked as "Estimation"). The dot-dashed lines incorporate the effect of "smooth 
accretion" (see text). Bottom panels: The same now for the stellar mass contribution from merging satellites. In each panel, the solid line 
shows the exact contribution from mergers of satellites into the final central galaxy, while the dashed line is the estimated contribution 
using the ZCZ07 approximation. 



contribution from smaller central galaxies in two ways: 
for a halo whose mass is around the peak of SFE or 
smaller, this assumption might overestimate the contri- 
bution, while for larger mass halos it can result in an 
underestimation (which will result in a stronger down- 
sizing behavior). 

The bottom panels of Figure [8] compare the rough esti- 
mation of the merger contribution from satellites to the 
stellar mass of the central galaxy (dashed lines) with the 
SAMs exact predictions (solid lines). The estimation is 
computed in the same halo mass range used in ZCZ07, 
to facilitate the comparison. In this case, we see that 



the approximation underestimates the merger contribu- 
tion of satellites to the final stellar mass, especially at 
the highest halo masses probed. This may arise due to 
the merging of additional satellites to the brightest one 
in each progenitor halo, and may again lead to a stronger 
downsizing effect. We note that the behavior of the esti- 
mated satellite merger contribution goes in the opposite 
direction than that of smaller central galaxies. This leads 
to a partial cancellation such that the estimate of the to- 
tal merger contribution is reasonable, which also implies 
that the inferred contribution from recent star formation 
is not strongly affected by the approximate nature of the 
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above method. 

It is hard to draw definitive conclusions from all these 
tests given the intrinsic uncertainties in both the SAMs 
and the ZCZ07 estimations. We clarify that the latter 
assumptions investigated here have to do with trans- 
forming the ZCZ07 measurement of the growth of stel- 
lar mass in central galaxies as a function of halo mass 
(their Fig. 8), which is robust, to the overall contribu- 
tion of mergers and star formation to the stellar mass 
assembly (their Fig. 9). More care should certainly be 
given to these assumptions, taking into account smooth 
accretion and incorporating better treatment of central 
and sate llite dynamics deter mined from analytic mod- 
els (e.g., iZentner et al.l I2005T) or from simulations (e.g . , 
IWang et al J 1200ft IWhite eTaLl [20071 IWake et al.l ;2008) . 
Nonetheless, we expect the qualitative results of ZCZ07 
to still be valid (or even somewhat strengthened with 
these corrections, as discussed above), and believe that 
such phenomenological methods can provide powerful 
constraints on theories of galaxy formation and evolu- 
tion. 

5. SUMMARY AND DISCUSSION 

In this paper we study theoretical predictions for 
the evolution of stellar mass in galaxies as a func- 
tion of their host halo mass, using semi-analytic galaxy 
formation models based on the Millennium simula- 
tion. We utilize two different SAM implementations , 
the MPA (iCroton et al.ll2006t iDe Lucia fe Blaizoll2007| ) 
and Durham (jBower et al.l 120061 ) models. We investi- 
gate the different contributions to the growth of stel- 
lar mass, and the role of mergers and star forma- 
tion in the stellar mass assembly from z ~ 1 to 
the present. Such an investigation with SAMs is 
timely and important as several recent studies have 
started to explore the galaxy-halo connection and re- 

lated inferences on galaxy evolutio n (e. g., ZCZ07 ; 

White et alJ[2007tlConrov fc Wechslerl[200l IWake et al 
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Leauthaud et al.l l201lT iNeistein et alJ l2011b| ). These 
studies employ different phenomenological approaches 
utilizing observed statistical properties of galaxies, such 
as correlation functions, abundances of galaxies and stel- 
lar mass functions. In our study, we particularly compare 
the SAM predictions to the methodology and results pre- 
sented in ZCZ07, to assess the potential of such studies 
to constrain galaxy formation models and to guide future 
efforts of modeling galaxy evolution. 

We find that the SFE, the fraction of baryon mass 
that has converted into stars in the central galaxies, as 
a function of halo mass has a peaked distribution, with 
a maximal value of ~ 23% at z ~ and ~ 18% at z ~ 
1. The location of the peak shifts toward lower halo 
mass with time, reflecting "halo downsizing" . Both SAM 
models produce similar results for the growth of stellar 
mass in central galaxies from z ~ 1 to as a function of 
the present-day halo mass. At both redshifts, the central 
galaxy stellar mass increases rapidly with halo mass for 
relatively low-mass halos (below - 2 x 10 12 /i _1 M© ) and 
at a lower rate for more massive halos. The fraction of 
stellar mass already in place at z ~ 1 also varies with final 
halo mass: it is about 50% (40%) for the MPA (Durham) 
model at the low-mass end, increasing with halo mass to 



about 65% for halos with mass ~ a few x10 12 /i _1 Mq, 
and decreases somewhat at the highest halo mass probed. 

The SAM predictions for the different contributions to 
the stellar mass assembly since z ~ 1 indicate that star 
formation is more important in low mass halos (~ a few 
xlO n /i _1 M0), while accretion through mergers domi- 
nates at the high-mass end (~ 10 13 /i _1 M Q ) , where star 
formation is negligible. In the intermediate regime both 
these processes contribute. This trend with halo mass is 
another manifestation of downsizing. The two SAMs pro- 
vide similar results, differing mostly in their predictions 
for the contribution of smaller central galaxies merging 
with the main central galaxy. This likely arises from dif- 
ferences in the galaxy formation prescriptions and in the 
merger trees of these models. 

We also study the predictions of the MPA SAM for the 
stellar mass growth since z ~ 2. The trends found are 
very similar to those for z ~ 1, including the presence of 
the "downsizing" pattern. As expected, much less stellar 
mass is already in place in the main progenitors com- 
pared to z ~ 1 and the contribution from merging of 
smaller central galaxies is considerably larger. Further- 
more, the contribution from star formation is important 
at all halo masses, even at the high-mass end. 

Our study is motivated by ZCZ07 who develop a novel 
phenomenological approach to study galaxy evolution by 
connecting galaxy clustering results at different epochs 
through the growth of the hosting halo mass. Such ap- 
plications can potentially provide important constraints 
for galaxy formation models as a function of the host halo 
mass, which is a fundamental parameter in such models. 
We compare our finding to those of ZCZ07. We find that 
the SAMs and ZCZ07 produce similar trends for the stel- 
lar mass assembly in halos, however, there are significant 
quantitative differences. The SFE of central galaxies as 
a function of halo mass at both epochs in the SAMs and 
ZCZ07 are qualitatively similar, with the same overall 
peaked shape and halo downsizing. The main discrep- 
ancy appears at z ~ 1 where the MPA SAM predicts a 
- 50% higher SFE than ZCZ07. 

The differences are also apparent when contrasting, for 
these two approaches, the stellar mass content of ha- 
los at the two epochs as a function of present-day halo 
mass (Fig. 5). While the overall trends are in qualita- 
tive agreement, there are striking differences. Again, the 
most significant difference is that the SAMs predict a 
larger stellar mass content at z ~ 1 for all halo masses 
(with a bigger discrepancy for lower- mass halos). The 
results from the SAMs and ZCZ07 are in better agree- 
ment at z ~ 0, however, the SAMs underestimate the 
stellar mass in central galaxies in present-day halos more 
massive than ~ 10 /i -1 Mq. Note that for halos larger 
than this characteristic mass, AG N feedback starts t o 
playing an important role (see e.g.. lStringer et aDl2008f l. 
suggesting its effect might be overestimated within these 
models, resulting in over-quenching stellar mass growth. 

When looking at the fraction of stellar mass in present- 
day central galaxies that is already in place at z ~ 1, 
there is a factor of two disagreement between the SAMs 
and the ZCZ07 results. For example, ZCZ07 obtain that 
about 30% of the stars in halos of ~ 3 x 10 12 /i~ 1 M Q are 
formed by z ~ 1, while the SAMs predict about 60%. 
The discrepancy is significant even when conservatively 
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accounting for a possible underestimation of the DEEP2 
stellar masses. Note also that these differences are al- 
ready present at higher redshifts, as we find that the 
ratio of stellar mass in central galaxies at z ~ 2 in the 
MPA SAM is comparable to the one predicted by ZCZ07 
at z r~j 1. Our results are in agreement with previous 
works that studied SAM predictions, whic h found an ex- 
cess of stars already in place by z ~ 1 (e.g.. lCroton et all 
2006; [Kitzbichler fc Whiti l2007t) . Those, however, were 
focused on integrated properties and not explicitly as a 
function of halo mass as we show here. 

The SAM predictions for the total amount of stars ac- 
quired through merging on top of that already in place 
at z ~ 1 are, at first order, in good agreement with 
ZCZ07 results. However, the individual contributions to 
the central galaxies stellar mass from mergers of smaller 
centrals and satellite mergers are markedly different. In 
ZCZ07 the significant merger contribution arises from the 
smaller central progenitors. In contrast, the SAMs pre- 
dict a large contribution from mergers of satellites. It is 
the partial cancellation of these opposing differences that 
leads to a reasonable agreement of the total merger con- 
tribution. As a whole, the SAMs and the ZCZ07 results 
lead to a similar behavior of the star formation contribu- 
tion with halo mass. 

While the comparison between the ZCZ07 observa- 
tional results and the SAMs predictions is informative, 
there are some simplified assumptions in the former. 
ZCZ07 apply the method as a proof of concept and point 
out that there is room for improvement with more sophis- 
ticated applications. With the SAMs, we are able to ex- 
amine different working assumptions employed by ZCZ07 
and guide these efforts. For instance, to link galaxies at 
two epochs, ZCZ07 use the average relationship between 
the present-day halo mass and the mass of the main pro- 
genitor at z ~ 1, neglecting the individual scatter among 
halos. We test the validity of this assumption, using the 
full assembly information available in the SAMs, finding 
that it results in negligible differences. 

On the other hand, some of the assumptions made by 
ZCZ07 to estimate the overall contribution from mergers 
and star formation can certainly be improved. In par- 
ticular, the original ZCZ07 estimation does not take into 
account the smooth accretion of dark matter particles to 
the final halo mass. In the Millennium Simulation the 
smooth accretion since z ~ 1 accounts for about ~ 30% 
of the final halo mass. It is difficult, however, to esti- 
mate the contribution of stellar mass from very small 
halos, below the resolution of current numerical simu- 
lations. If smooth accretion is as significant in the real 
universe, it will certainly be needed to be included in such 
approximations. More realistic SFE dependence on halo 
mass can also be implemented to improve the estimation 
method. 

Conroy & Wechsler (2009) present related calcula- 
tions for the evolution of stellar masses and star forma- 
tion. They use abundance matching to monotonically 
link galaxies to halos. They predict similar, but more 
dramatic trends than both the SAMs and the ZCZ07 
approach. For instance, they suggest that galaxies in 
lower mass halos (~ 10 11 /i _1 M Q ) grow their stellar mass 
purely by star formation, while essentially all of the mass 



is already in place by z ~ 1 in present-day halos of 
lO 13 /i _1 M0. Their results support a picture in which 
stellar mass grows only via star formation, suggesting 
that the stellar mass from smaller progenitors does not 
merge into the central galaxy, but remains as satellites 
or diffuse light. Additional studies are needed to fully 
clarify their differences with the results prese nted here 
and in ZCZ07. Recently, iNeistein etaLl (|2011al) has crit- 
ically studied the assumptions made in the abundance 
matching method using SAM catalogs and found impor- 
tant differences, indicating that environmental processes 
may be important. 

We have demonstrated that phenomenological meth- 
ods such as ZCZ07 are powerful for studying galaxy for- 
mation and evolution. They provide key constraints for 
theoretical models, such as the SAMs, as a function of 
halo mass. By highlighting remaining shortcoming of 
galaxy formation models, they can guide to improving 
theoretical predictions at high redshift and increase our 
understanding of the complex picture of galaxy formation 
and evolution. At the same time, while ZCZ07 is useful 
as a proof of concept, future work should use more sophis- 
ticated methods applied to better data, and the SAMs 
can serve an important role in developing and testing 
such methods. 

Future work will also explore the role of environment 
in the buildup of stellar mass within the host halos. One 
of the main assumptions in the current HOD frame- 
work is that the galaxy content in halos depends only 
on the halo mass, and is independent of the large-scale 
environment where the halo is located. Recent theoreti- 
cal studies have shown that the clustering properties of 
dark matter halos depend on the large-scale environment 
(the s o-called halo ass embly bias;lGao. Springel fc White! 
1 20051: iWechsler et all 120061 iCroton. Gao fc White! 120071 : 
Ping. Suto fc Mdl2007l ). This environmental effect might 
also impact galaxy properties and galaxy clustering 
(|Zhu et al.|[200l IZu et at] [2008). Using the SAMs, we 
may be able to test the effect of large-scale environ- 
ment on stellar mass assembly and galaxy evolution (c.f., 
IHovle. Jimenez fc Verdd[2011h . Moreover, we could use 
SAM results to incorporate environmental effects to phe- 
nomenological methods such as ZCZ07, increasing the 
constraining power of galaxy clustering data on galaxy 
formation models. 
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